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ABSTRACT 

We analyze high-resolution spectra of the nearby (1.34 pc) stars a Cen A (G2 V) and a Cen B (K1 V), 
which were obtained with the Goddard High Resolution Spectrograph on the Hubble Space Telescope . 
The observations consist of echelle spectra of the Mg n 2800 A and Fe n 2599 A resonance lines and the 
Lyman-a lines of hydrogen and deuterium. The interstellar gas has a velocity (r = —18.0 + 0.2 km s" *) 
consistent with the local flow vector proposed for this line of sight by Lallement & Bertin (1992). The 
temperature and nonthermal velocity inferred from the Fe ii, Mg ii, and D i line profiles are 
7 = 5400 + 500 K and £ = 1.20 ± 0.25 km s \ respectively. However, single-component fits to the H i 
Lyman-a lines yield a Doppler parameter (h Hl = 11.80 km s' 1 ) that implies a significantly warmer tem- 
perature of 8350 K, and the velocity of the H i absorption (i? = —15.8 ±0.2 km s -1 ) is redshifted by 
about 2.2 km s" 1 with respect to the Fe ii, Mg ii, and D i lines. The one-component model of the 
interstellar gas suggests log N H[ = 18.03 ± 0.01 and D/H = (5.7 ± 0.2) x 10 6 . These parameters lead to 
a good fit to the observed spectra, but this model does not explain the higher temperature and redshift 
of H i relative to the other interstellar lines. 

The most sensible way to resolve the discrepancy between H i and the other lines is to add a second 
absorption component to the H i lines. This second component is hotter ( T ^ 30,000 K), is redshifted 
relative to the primary component by 2-4 km s _1 , and has a column density too low to be detected in 
the Fe ii. Mg ii, and D i lines. We propose that the gas responsible for this component is located near 
the heliopause, consisting of the heated H i gas from the interstellar medium that is compressed by the 
solar wind. This so-called “hydrogen wall” is predicted by recent multifluid gasdynamical models of the 
interstellar gas and solar wind interaction. Our data provide the first measurements of the temperature 
and column density of H i in the hydrogen wall. After considering the effects that a corresponding 
hydrogen wall around a Cen would have on our analysis, our best estimates for the parameters of the 
solar hydrogen wall are log = 14.74 ± 0.24, b { ^\ = 21.9+1.7 km s 1 (corresponding to 7 = 
29,000 ± 5000 K), and v\}\ > - 16 km s" *. 

Unfortunately, the existence of this heated H i reduces our ability to compute the H i column density 
of the interstellar medium accurately because, with slight alterations to our assumed stellar Lyman-a 
profiles, we discovered that acceptable two-component fits also exist with log N Hl ~ 17.6. We, therefore, 
quote large error bars for the H i column density along the a Cen line of sight, log N HI = 17.80 ± 0.30. 
For this range in iV HI , n H j = 0.15 cm -3 (± a factor of 2) and D/H = (0.5- 1.9) x 10 5 . This is the first 
direct measurement of the H i density in a local cloud and allows us to predict the distance from the 
Sun to the edge of the local cloud along various lines of sight. This range in D/H is consistent with the 
value D/H = 1.6 x 10' 5 previously derived for the Capella and Procyon lines of sight. We cannot tell 
whether the D/H ratio varies or is constant in the local interstellar medium, but we do find that the 
D i/Mg n ratio for the a Cen line of sight is about 4 times smaller than for the Capella and Procyon 
lines of sight. Therefore, either D/H or the Mg depletion varies significantly over distance scales of only 
a few parsecs. 

Subject headings : ISM: abundances — ISM: kinematics and dynamics — radio lines: ISM 
stars: individual (a Centauri) — ultraviolet: ISM 


L INTRODUCTION 

We analyze in this paper our Goddard High Resolution 
Spectrograph (GHRS) observations of interstellar gas 
absorption along the line of sight toward both components 
of the a Cen visual binary system — a Cen A (G2 V) and a 
Cen B (K1 V). Our objective was to measure very accurately 
the cosmologically important deuterium/hydrogen (D/H) 


1 Based on observations with the NASA/ESA Hubble Space Telescope , 
obtained at the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. 

2 JILA, University of Colorado, Boulder, CO 80309-0440. 


ratio in a very short (1.34 pc) line of sight that we antici- 
pated should have rather simple and homogeneous physical 
properties. This program is one component of a larger study 
of D/H in the local interstellar medium (LISM). Our pre- 
vious studies of the 12.5 pc line of sight toward Capella 
(Linsky et al. 1993, hereafter Paper I) and the 3.5 pc line of 
sight toward Procyon (Linsky et al. 1995, hereafter Paper 
II) led us to conclude that D/H = (1.6 + 0.09 [and +0.05, 
— 0.10 systematic error]) x 10' 5 in the LISM. Because 
deuterium is converted by nuclear reactions in stars to 3 He 
and 4 He (a process usually called “astration ”) over the life- 
time of the Galaxy, and this deuterium-depleted gas is 
returned to the interstellar medium (ISM) by stellar winds 
and supernova explosions, the D/H ratio in the LISM pro- 
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vides a lower limit to the cosmologically important primor- 
dial D/H ratio and also provides an important constraint 
on Galactic chemical evolution models. By studying the gas 
in the a Cen line of sight, we are testing whether the D/H 
ratio is essentially constant in the LISM, as implied by the 
short timescale for mixing of interstellar gas, or whether 
there are significant variations in the D/H ratio, as sug- 
gested by Copernicus observations (see, e.g., Dupree, 
Baliunas, & Shipman 1977; Bruston et al. 1981). 

The high spectral resolution (A/AA = 90,000) and signal- 
to-noise (S/N) capabilities of the GHRS and the suppress- 
ion of geocoronal Lyman-a emission achievable with the 
use of the 0725 x 0725 Small Science Aperture (SSA) make 
the GHRS a far superior instrument than Copernicus or the 
International Ultraviolet Explorer (WE) for studies of D/H 
and the dynamics of the LISM. The GHRS has been used 
so far to study a number of lines of sight through the LISM 
including as targets Capella, Procyon, G191— B2B 
(Lcmoine et al. 1995), Sirius (Lallement et al. 1994), a PsA 
(Ferlet et al. 1995), Altair, Vega, /? Pic, 8 Cas (Lallement et 
al. 1995), and € CMa (Gry et al. 1995). 

Analyses of these lines of sight have led to the following 
picture of the LISM as summarized by Lallement et al. 
(1995). The Sun lies near the edge of what is called the local 
interstellar cloud (LIC), which is also called the AG cloud 
because it is centered in the anti-Galactic center direction. 
This cloud appears to move as a rigid body because high- 
resolution spectra of nearly all nearby stars show significant 
interstellar absorption at velocities within 1 km s _1 of the 
projected velocity of the LIC How vector. In addition, most 
stars show additional absorption features centered at other 
velocities, which is usually interpreted as evidence of other 
clouds that presumably lie outside the LIC. The gas tem- 
perature and nonthermal motions (often called turbulence) 
in the LIC can be derived from the line widths of ions that 
differ greatly in atomic mass. The resonance lines of D i, 
Fe it, and Mg n arc quite suitable for this purpose. The 
Capella and Procyon lines of sight (see Paper II) have pro- 
vided perhaps the most accurate values for these param- 
eters, T = 7000 + 500 (and ±400 systematic error) K and 
£ = 1.20 ± 0.27 km s’ 1 , where £ is the most probable speed 
for the assumed Gaussian nonthermal motions. Frisch 
(1995) provides a comprehensive review of the physical and 
dynamic properties of the LISM. 

The a Cen line of sight provides a critical test of this 
picture because it is the only star studied so far by the 
GHRS that lies in the general direction of the Galactic 
center, and the interstellar Fe n and Mg ii absorptions for 
this line of sight have velocities inconsistent with the LIC 
flow vector (Lallement et al. 1995). In fact, Lallement & 
Bertin (1992) placed a Cen in a different cloud, their so- 
called G (for Galactic) cloud, with a slightly different flow 
vector. Analysis of the a Cen spectra may, therefore, provide 
answers to the questions of w here the LIC ends and the G 
cloud begins in this direction, the physical properties of the 
G cloud, and whether D/H in the G cloud differs from that 
in the LIC. The GHRS has measured D/'H so far only for 
gas in the LIC. Most stars studied to date with the GHRS 
probably lie outside the LIC, and thus the derived H i 
column densities provide only a lower limit to the H i 
number density in the LIC. On the other hand, the line of 
sight to a Cen likely lies entirely or mostly within the G 
cloud, so that one can determine the H i number density in 
the G cloud from the H i column density. 


When analyzing lines of sight with low H i column den- 
sities, one should consider whether hydrogen in the helio- 
sphere contributes significantly to the total hydrogen 
absorption. The interaction between the ionized solar wind 
and the partially ionized interstellar gas flowing toward the 
Sun is now thought to create three discontinuities in the 
flow: a bow shock located at 200 500 AU in the upwind 
direction, the heliopause separating the flow of interstellar 
plasma from that of the solar-wind plasma, and the solar- 
wind termination shock located at about 100 200 AU in 
the upwind direction where the solar wind makes the tran- 
sition from supersonic to subsonic flow (cf. Baranov & 
Malama 1995). The neutral component of the interstellar 
gas flows through the solar system, but charge exchange 
with solar-wind protons slows and heats this flow, forming 
a region of enhanced neutral-hydrogen density called the 
“hydrogen wall” located between the bow shock and the 
heliopause at roughly 200 AU in the upwind direction. We 
believe that we have detected the hydrogen wall through 
redshifted absorption by the heated hydrogen in the a Cen 
line of sight, and for the first time we can estimate its 
column density and temperature (see § 4.5). 

The determination of D/H in the LISM involves the mea- 
surement of the neutral-deuterium column density (N Dl ) 
and the neutral-hydrogen column density ( N H] ) by inter- 
stellar absorption in the hydrogen and deuterium Lyman-a 
lines against the stellar chromospheric Lyman-a emission 
line. Since the D i Lyman-a line is optically thin toward a 
Cen and relatively narrow, we can interpolate the stellar 
emission line background to obtain an accurate value of 
iV D1 . On the other hand, the H l Lyman-a absorption 
feature is very broad and saturated, so we must model the 
unknown stellar emission line. Errors in its assumed line 
shape can translate into large systematic errors in the 
derived value of N Hl . Thus, the main uncertainty in D/H is 
the H i column rather than the D i column. For this reason, 
we pay very careful attention to the roles played by the 
uncertain stellar emission line and the influence of helio- 
spheric hydrogen. This care required a great many simula- 
tions to characterize the random and systematic errors to be 
placed on the derived parameters. In this paper, we consider 
parameter errors to be random when they characterize the 
quality of the fit to the data and errors in the derived 
parameters to be systematic when they arc a consequence of 
the unknown emission-line profile background. 

The GHRS spectra are described in § 2, and we present in 
§ 3 our determination of the interstellar parameters and 
D/H ratio under several different assumptions. In § 4, we 
compare our derived interstellar parameters with those 
obtained for other lines of sight, and we also discuss models 
of H i in the heliosphere and how they relate to our analysis. 

2. GODDARD HIGH RESOLUTION SPECTROGRAPH 
OBSERVATIONS OF ALPHA CENTAURI 

We observed a Cen A on 1995 May 1 and a Cen B on 
1995 May 5 by using the echelle gratings of the GHRS and 
the SSA to obtain the highest possible spectral resolution 
and to suppress the geocoronal Lyman-a emission. For a 
description of the GHRS instrument, see Brandt et al. (1994) 
and Heap et al. (1995). Soderblom et al. (1994) describe the 
improved characteristics of the GHRS with COSTAR when 
a Cen was observed. Table 1 summarizes the observations, 
including the spectral resolutions and integration times. 
The echelle data were obtained in the FP-SPLIT mode to 
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TABLE 1 

Summary of Goddard High Resolution Spectrograph Observations 


Spectral Spectral Exposure Start 

Aperture and Range Resolution Time Time 

Grating Substep Pattern (A) (km s _1 ) (s) (UT) Important Features 


a Cen A (1995 May 1) 

G140L LSA 5 1200-1490 360 186 14:47 H i 1216 A 

EA-46 SSA 9 1212 1219 3.57 2786 15:02 Hi + Di 1216 A 

EB-20 SSA 7 2792-2806 3.54 485 17:40 Mg ii 2796, 2803 A 

EB-22 SSA 7 2593-2605 3.27 646 17:52 Fe ii 2599 A 


a Cen B (1995 May 5) 

EB-20 SSA 7 2792-2806 3.54 700 10:25 Mg ii 2796, 2803 A 

EA-46 SSA 9 1212-1219 3.57 3231 12:04 H I + D i 1216 A 


reduce fixed-pattern noise and were processed with the 
version of the CALHRS software available in 1995 May. 
The individual readouts of the FP-SPLIT spectra were 
combined by using a cross-correlation procedure called 
HRS_MERGE (Robinson et al. 1992). Images of the plati- 
num calibration lamp were taken prior to all observations 
and used to calibrate the wavelengths of the spectra. 

We observed the Mg ii h and k resonance lines (at 2803 A 
and 2796 A) in both stars and the Fe ii multiplet UV1 (at 
2599 A) in a Cen A only. Our objective in observing lines of 
these elements, which have much larger atomic weights 
than H or D and therefore have narrower absorption lines, 
was to search for multiple velocity components in this line 
of sight and to measure the thermal and turbulent broaden- 
ing independently. Figure 1 shows the complete echelle B 
spectra. The high resolution of the Ech-B grating separates 
the very narrow and deep interstellar Mg ii and Fe ii lines 
from the broader stellar chromospheric Mg n emission lines 
and photospheric Fe ii absorption line (see § 3.1). 

Our spectra of the Mg ii h and Fe ii lines of a Cen A are 
essentially identical in shape to the GHRS spectra obtained 
in 1992 by Lallement et al. (1995), except that our fluxes are 


lower (by factors of 5.6 and 4.8 for the Mg ii and Fe ii 
spectra, respectively), presumably because of the star’s being 
misplaced at the edge of the small aperture. Comparison 
with Mg ii fluxes measured by IUE (see, e g., Cerruti-Sola, 
Cheng, & Pallavicini 1992) confirms that the fluxes of the 
1992 observations are correct, so we have normalized the 
fluxes of our spectra to agree with those of the 1992 spectra. 
(Fig. 1 and all subsequent figures show the corrected fluxes.) 
Acquisition and precise centering in the 0"25 x 0'.'25 SSA 
for the a Cen stars are difficult because of the rapid motion 
of these stars, the brightness of each star, and the proximity 
of the bright companion star (about 18" away). Our wave- 
length calibrations are accurate because our a Cen A and a 
Cen B spectra and the 1992 a Cen A spectra show the same 
ISM velocities for the Fe ii and Mg ii lines. 

Figure 2 shows the spectra of the Lyman-a region of both 
stars. Superposed on the broad Lyman-a emission lines 
formed in the chromospheres of the stars are broad saturat- 
ed absorption features (centered at 1215.6 A) due to inter- 
stellar H i and narrower absorption lines due to interstellar 
D i, predicted to lie at —0.3307 A (corresponding to —81.6 
km s _1 ) relative to H i. The flux of the Lyman-a line of a 
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Fig. 1. — Our echelle-B observations of the Mg ii h and k lines (at 2802.7 and 2795.5 A) for a Cen A and a Cen B ( upper two panels) and the Fe n 2599.4 A 
line of a Cen A (bottom panel). The dotted lines mark the locations of interstellar absorption lines. 
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Fig. 2. — Our echelle-A observations of the Lyman-a lines of a Cen A 
and a Cen B, showing broad interstellar H i absorption centered at 1215.6 
A and narrow D i absorption near 1215.25 A. 


Cen A was a factor of 3.3 times lower than that observed 
during the G140L observation listed in Table 1. This con- 
firms that a Cen A was not well centered in the SSA during 
our high-resolution observations. Since the G140L spec- 
trum was taken through the Large Science Aperture (LSA) 
rather than the SSA, the flux calibration should be much 
more reliable. Therefore, we increased the fluxes of our 
Lyman-a spectrum of a Cen A by a factor of 3.3 to obtain 
consistency with the G140L data. Since the velocity of the 
interstellar D i line is consistent with the velocities of the 
Fe ii and Mg ii lines, we believe the wavelength calibration 
is accurate. 

Owing to a guide star acquisition failure, the a Cen B 
observations were obtained in coarse track guiding mode 
on a single guide star, meaning the accuracy of the centering 
of the star in the aperture is highly suspect for these obser- 
vations. The first science observation after acquisition was 
the Mg ii observation. We find that the Mg ii fluxes we 
observed are consistent with those observed by the IUE 
(see, e.g., Cerruti-Sola et al. 1992). We also find that the 
velocities of the interstellar Mg ii absorption lines are the 
same as those seen in the a Cen A spectra, which means that 
the wavelength calibration of this spectrum is accurate. 
Because the fluxes and wavelengths appear to be correct, a 
Cen B was apparently well centered in the SSA at the begin- 
ning of our observations. However, during the subsequent 
observation of the Lyman-a line we believe the star must 
have been slowly drifting out of the aperture because 
inspection of the 16 individual readouts of this FP-SPLIT 
observation shows that the count rates decreased through- 
out the observation. We therefore use the first readout of 


the observation to flux-calibrate the spectrum. This pro- 
cedure amounts to multiplying the fluxes of the combined 
spectrum by a factor of 1.3. We also find that the velocity of 
the interstellar D i line is redshifted by 0.9 km s _1 with 
respect to the interstellar D i absorption line observed for a 
Cen A (which has a velocity consistent with the velocities 
measured for Fe ii and Mg ii). This discrepancy is also 
probably due to the star’s drifting out of the aperture during 
the observation, and so we have changed the wavelengths 
by 0.9 km s - 1 to make this observation consistent with the 
others. (Note that Fig. 2 and all subsequent figures show the 
corrected fluxes and wavelengths for the Lyman-a lines of a 
Cen A and a Cen B.) 

The D i lines in Figure 2 are well separated from the 
interstellar H i absorption, and since they are not very deep 
they must be optically thin. On the other hand, the inter- 
stellar H i lines are very saturated with sharp core bound- 
aries, indicating that these lines lie on the flat part of the 
curve of growth with H i column densities log /V H1 ~~ 18. 
Thus, the absorption line shape depends both on the H i 
column and the broadening parameter b Hl . The saturated 
H I absorption core has positive flux, which can result only 
from scattered light in the GHRS spectrograph. This 
problem has been seen previously in our Lyman-a spectra 
of Capella (Paper I). In our analysis and in all figures except 
Figure 2, we correct for this scattered light by simply sub- 
tracting the mean flux level in the saturated core 
(9.9 x 10“ 13 erg cm - 2 s“ 1 A - 1 and 1.9 x 10" 12 erg cm 2 
s _1 A -1 for a Cen A and a Cen B, respectively) from the 
fluxes observed between 1215.35 and 121 5.85 A. 

3. PROPERTIES OF THE LINE OF SIGHT TO 
ALPHA CENTAURI 

3.1. Single-Component Fits to the Line Profiles 

We first analyze each line profile of each star separately 
without imposing any constraints on the derived param- 
eters, except to minimize the x Z of the fit. It is important, 
however, to estimate as accurately as feasible the intrinsic 
stellar emission or absorption line upon which the inter- 
stellar absorption is superposed, as uncertainties in the 
assumed stellar line profile can be the dominant source of 
error, especially for the highly saturated H i line. 

Figure 3 shows our fit to the a Cen A stellar and inter- 
stellar Fe ii 2599 A lines, for both our data and the higher 
S/N archival data. Figures 4 and 5 show our fits to the Mg n 
h and k lines for both stars and, for the a Cen A h line, both 
data sets. The thin solid lines in the figures represent our 
best guesses for the intrinsic stellar line profiles, estimated 
by using polynomial fits to the spectral regions blueward 
and redward of the interstellar absorption lines. The thick 
solid lines, which fit the observations very well, are the con- 
volution of the computed interstellar absorption lines 
(dotted lines) and the instrumental point-spread function, 
which is approximated as a Gaussian with FWHM — 3.7 
pixels (Gilliland 1994). There is no evidence for any asym- 
metry in the absorption lines, and the residuals of our one- 
component fits (shown below each spectrum) reveal no 
dependence on velocity that would suggest the presence of a 
second component. 

The parameters for these fits are listed in Table 2 together 
with random errors estimated by using Monte Carlo simu- 
lations. We have matched the observed profiles very well, 
but the derived ISM parameters for the two Mg n lines and 
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Fig. 3. — Our best fits to the Fe n interstellar absorption lines observed 
toward a Cen A in 1992 and in 1995. The parameters of the fits are given in 
Table 2. The data are shown in histogram form. The thin solid lines are our 
estimates for the stellar emission above the absorption lines, computed by 
using polynomial fits to spectral regions blueward and redward of the 
absorption. The dotted line is the best-fit absorption line before convolu- 
tion with the instrumental profile, and the thick solid line is this absorption 
line after instrumental broadening. The residuals are plotted below each fit. 

for the same lines in the two stars often lie outside each 
other’s formal statistical errors. This is probably due to the 
systematic uncertainties in our estimates of the intrinsic 
stellar line profiles in the vicinity of the interstellar features. 
The fr FeU -values for the two data sets differ from each other 
by a little more than 2 tr, but the Fe ii lines have widths very 
close to the instrumental resolution, making it difficult to 
infer an accurate value for b FeU , especially for the noisier 
1995 data. The heliocentric velocities listed in Table 2 are 
consistent with the predicted — 18.1 km s“ 1 velocity of the 
local interstellar flow for the Galactic center direction pro- 
posed by Lallement & Bertin (1992). 
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Fig. 4. — Our best fits to the Mg ii interstellar absorption lines observed 
toward a Cen A in 1992 and in 1995. The parameters of the fits are given in 
Table 2. The various lines that are displayed have the same meaning as in 
Fig. 3. 


In Figure 6, we fit the interstellar D i absorption lines by 
using polynomial fits to estimate the “continua” against 
which the D i absorption lines are measured. Each of these 
continua actually represents the combination of the stellar 
emission line plus some absorption in the wings of the inter- 


TABLE 2 


Interstellar Parameters for the Single-Component Fits 


Star 

Ion 

Line 

Velocity 
(km s 1 ) 

b 

(km s~ l ) 

log N 

X 

x 2 

a Cen A a 

Fe ii 

2599.396 

-18.2 + 0.1 

1.43 + 0.02 

12.441 ±0.004 

1.64 

0.701 

a Cen A 

Fe il 

2599.396 

-17.7 + 0.1 

1.78 + 0.14 

12.455 + 0.022 

1.39 

1.641 

a Cen A 

Mg ii 

2795.528 

-18.1 + 0.1 

2.32 + 0.04 

12.698 + 0.017 

5.50 

1.708 

a Cen A a 

Mg ii 

2802.705 

-17.8 ± 0.1 

2.26 + 0.02 

12.691 +0.003 

2.76 

5.000 

a Cen A 

Mg ii 

2802.705 

-17.8 T 0.1 

2.21 + 0.04 

12.776 + 0.010 

3.35 

1.624 

a Cen B 

Mg ii 

2795.528 

-18.2 + 0.1 

2.49 + 0.01 

12.640 + 0.003 

4.48 

5.349 

a Cen B 

Mg ii 

2802.705 

-18.0 + 0.1 

2.19 + 0.01 

12.802 + 0.003 

3.61 

4.924 

a Cen A 

D i 

1215.339 

-18.2 + 0.1 

6.55 + 0.27 

12.799 + 0.015 

0.72 

0.928 

a Cen B 

D i 

1215.339 

-18.2 + 0.1 

6.88 + 0.17 

12.775 + 0.009 

0.65 

1.334 

a Cen A 

H l 

1215.670 

-15.9 ± 0.1 

11.80 + 0.05 

18.031 +0.003 

68671 

1.090 

a Cen B 

H i 

1215.670 

-15.6 + 0.1 

11.81 ±0.02 

18.027 + 0.002 

68111 

2.211 


Spectra obtained by Lallement et al. 1995. 
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Fig. 5. — Our best fits to the Mg n interstellar absorption lines observed 
toward a Cen B in 1995. The parameters of the fits are given in Table 2. 
The various lines that are displayed have the same meaning as in Fig. 3. 



Velocity (km s' 1 ) 

Fig. 6 . — Our best fits to the D i interstellar absorption lines observed 
toward a Cen A and a Cen B. The parameters of the fits are given in Table 
2. The various lines that are displayed have the same meaning as in Fig. 3. 
The dotted lines showing the best-fit profile before instrumental smoothing 
are barely distinguishable from the thick solid lines. 


stellar H i absorption profile. For Capella and Procyon it 
was not possible to fit D i separately from H i (see Papers I 
and II), as we have done here, because for each of those two 
cases the widths of the Lyman-a emission line and the inter- 
stellar H i absorption line act to create a local flux 
maximum somewhere in the spectral region of the D i line. 
This local flux maximum makes it very difficult to estimate 
the “ continuum ” for the D i absorption line by looking 
only at that spectral region. For a Cen A and a Cen B, 
though, we believe that we can accurately extrapolate a 
continuum above the D I absorption lines without deter- 
mining the shape of the stellar Lyman-a profile and then 
correcting for the interstellar H i absorption. For all of our 
fits to the D i and H i lines, we have included both com- 
ponents of the fine structure doublet, which has a velocity 
spacing of 1 .33 km s " 1 . 

In fitting the H I line profiles, we initially ignore the D i 
lines by using the same polynomial interpolation scheme 
used for the D i fits to remove the interstellar D i absorp- 
tion. The process we use to construct the intrinsic stellar 
Lyman-a profiles is described below and illustrated in 
Figure 7. Because both a Cen A and a Cen B are solar-like 
main-sequence stars, we use the quiet-Sun profile from 
Brekke et al. (1991) as the starting point of our analysis 
(dotted lines in the figure). However, a Cen A and a Cen B 



Fig. 7. — An illustration of how we estimated the intrinsic stellar 
Lyman-a profiles of a Cen A and a Cen B (see text for details). The dotted 
tine is the quiet-Sun Lyman-a profile taken from Brekke et al. (1991). The 
dashed lines are our initial guesses for the stellar Lyman-a profiles. For a 
Cen A (a Cen B), the wings of the dashed line are based on broadening 
(narrowing) the solar Lyman-a profile by 10% (6%). The centers of the 
dashed lines are based on the Mg it k lines observed for a Cen A and a Cen 
B. The assumed stellar profiles were altered to improve the quality of 
single-component fits to the interstellar H i absorption lines, and the solid 
lines show the stellar profiles after these alterations. 
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have luminosities different from the Sun’s, and the widths of 
optically thick chromospheric lines such as Lyman-a are 
dependent on the stellar luminosity, with the higher lumi- 
nosity stars having the broader lines. This phenomenon, 
which is called the Wilson-Bappu effect (Wilson & Bappu 
1957), was first studied using the Ca n H and K lines, but a 
similar correlation exists for the Mg n h and k lines (see, e.g., 
Vladilo et al. 1987) and the Lyman-a line (Landsman & 
Simon 1993). We have compared the Mg ii k line profile of 
the Sun with those of a Cen A and a Cen B, and we find that 
the a Cen A line is roughly 10% broader than the solar line 
and the a Cen B line is roughly 6% narrower. This result is 
consistent with the Wilson-Bappu effect, since a Cen A is 
somewhat more luminous and a Cen B is less luminous than 
the Sun. Assuming that the widths of Lyman-a scale like 
those of Mg ii, we next broaden (narrow) the solar Lyman-a 
line by 10% (6%) and scale the resulting profile to agree 
with the far wings of the a Cen A (a Cen B) line. This 
procedure provides our initial estimates for the wings and 
sides of the Lyman-a lines. 

The shapes of the centers of the Lyman-a profiles are 
estimated from the observed Mg ii k lines, because for the 
Sun and Procyon (see Paper II) we have found that the 
Mg ii k and Lyman-a lines have very nearly the same shape 
when plotted on a A2 scale (as opposed to a velocity scale). 
This similarity in shape can be understood in terms of both 
lines being optically thick chromospheric lines with self- 
reversals. We therefore use the observed a Cen A and a Cen 
B Mg ii k line profiles to estimate the shapes of the cores of 
the stellar Lyman-a lines, which we combine with the line 
wings estimated from the scaled solar Lyman-a profiles. The 
dashed lines in Figure 7 show the resulting profiles, which 
represent our initial guesses for the a Cen A and a Cen B 
Lyman-a profiles. Our initial fits to the data using these 
profiles were not very good, so we used the residuals of 
those first fits to alter the assumed profiles slightly and then 
we refit the data. We went through this procedure one more 
time before converging on the profiles shown as solid lines 
in Figure 7. The best fits to the data assuming these profiles 
are shown in Figure 8, and the parameters of the fits are 
given in Table 2. We note that the Lyman-a fluxes of the 
stellar Lyman-a profiles shown in Figures 7 and 8 are very 
similar to those that would be observed from the Sun at the 
distance of a Cen. The intrinsic central reversals are located 
entirely within the saturated interstellar absorption and 
thus play no role in the derived interstellar parameters. 

The fits displayed in Figure 8 match the data well, 
although the residuals reveal systematic discrepancies seen 
for both stars (near 1215.42 A and 1215.75 A) that cannot be 
removed by any reasonable alteration of the assumed stellar 
Lyman-a lines. Note that the computed Lyman-a line fits 
both the steep sides of the absorption core, which depend 
on both N u , and b H ly and the far wings of the line beyond 
the emission peaks, which depend only on N HI . It is very 
encouraging that the parameters of our single-component 
fits to the a Cen A and a Cen B data are almost identical 
(see Table 2). This is strong evidence that extracting infor- 
mation on the ISM from the observed H i absorption lines 
is a well-posed problem with a unique solution— at least for 
single-component fits. 

The line width parameter b = (0.0165T/A 4- £ 2 ) 1/2 km 
s' 1 , where A is the atomic weight and £ is the most prob- 
able speed (in km s” 1 ) of the assumed Gaussian turbulent 
motions. The combination of the 6-values for D (A = 2), 
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Fig. 8. — Our best single-component fits to the H l absorption lines of a 
Cen A and a Cen B. The parameters of the fits are given in Table 2. The 
data are shown in histogram form, where the D i absorption at 1215.25 A 
has been removed by using polynomial fits. The thin solid lines are our best 
estimates of the intrinsic stellar line profiles (see text and Fig. 7). The thick 
solid lines are the fits to the data, and the residuals are shown below each 
fit. The difference between the computed absorption lines before and after 
smoothing with the instrumental profile is too small to be visible on these 
plots. 


Mg (A = 24.3), and Fe (A — 55.8) allows us to solve for T 
and £ separately. The shaded region of Figure 9 shows the 
allowed range of these parameters given the values and 
uncertainties of b DU 6 Mgn , an d ^Fe ii for the two stars. From 
Figure 9, we quote values of T = 5400 + 500 K and 
£ = 1.20 + 0.25 km s _1 for the temperature and non- 
thermal velocity, respectively, for the a Cen line of sight. 

There are three surprising results from our one- 
component analysis of the a Cen ISM absorption lines that 
concern us and will be addressed in the next sections, where 
we explore models with multiple components and different 
intrinsic stellar emission line profiles. 

1. The ISM velocities measured by the fits to the Fe ii, 
Mg ii, and D i lines are nearly identical and collectively 
imply v = — 18.0 + 0.2 for the a Cen line of sight, but the 
centroid of the H i absorption is shifted by about +2.2 km 
s' 1 for both a Cen A and a Cen B. One expects, perhaps 
naively, that hydrogen and deuterium should be well mixed 
in the ISM and share the same flow velocities within mea- 
surement errors. However, Landsman et al. (1984) had pre- 
viously called attention to a similar discrepancy in their 
analysis of less accurate Copernicus and IUE spectra. They 
found that for a Cen A, the D i line is blueshifted by 8 + 2 
km s “ 1 relative to H i (or H i is redshifted by the same 
amount relative to D i). 

2. Second, the value of b H , and the temperature it implies 
are larger than the values inferred from the other line widths 
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Fig. 9. — Combination of the ^-values given in Table 2 for D i, Mg n, and Fe n allows us to solve for the temperature of the interstellar medium, T, and the 
nonthermal velocities present in the ISM, The allowable ranges of T and £ are shown for D i, Mg n, and Fe n separately, and the shaded region shows 
where these regions overlap. Also shown are the obviously discrepant allowable ranges of T and £ suggested by the single-component fits to H i (see § 3.1). 


(see Fig. 9). By using the values T = 5400 K and £ = 1.2 km 
s“ 1 derived from the D i, Mg n, and Fe ii lines, we compute 
6 hi = 9.51 km s' 1 , which is 2.3 km s _I smaller than 
observed. An alternative way of looking at this result is to 
assume the measured values for b H , and £, = 1.20 km s~ 1 to 
obtain T = 8350 K, almost 3000 K hotter than inferred 
from the other lines. Landsman et al. (1984, 1986), in their 
analysis of the Copernicus and IUE spectra, had also called 
attention to the very broad hydrogen profiles for these stars 
with b H , > 1 1 km s ' 1 for a Cen A and b H , > 1 3.8 km s “ 1 
for a Cen B. 

3. From the fits in Figures 6 and 8, we find that the D/H 
ratio is (5.86 + 0.24) x 10" 6 for the a Cen A data and 
(5.60 + 0.14) x 10' 6 for a Cen B. These D/H ratios are 
consistent within their cited errors and together imply a 
value of D/H = (5.7 + 0.2) x 10' 6 . This result is very differ- 
ent from the D/H = 1.6 x 10' 5 value found for the Capella 
and Procyon lines of sight. 

Are these results real or spurious, perhaps produced by 
inaccurate intrinsic Lyman-a profiles and/or by the 
assumption that the gas along this short line of sight is 
homogeneous? 

If true, the D/H = (5.7 + 0.2) x 10" 6 result could have 
profound implications for both cosmology and our under- 
standing of the LISM because this is a very different result 
from that found in Papers I and II (D/H = 1.6 x 10' 5 ). It is 
important, therefore, to determine whether our a Cen data 
are really inconsistent with D/H = 1.6 x 10“ 5 . Since D i is 
optically thin, log N DI should be well determined, and so 
decreasing log JV HI is the only possibility. For D/H to be 
1.6 x 10" 5 , log N h , must be decreased to 17.6. In Figure 10 
we have fit the Lyman-a lines again, but this time we have 


forced log N Hl to be 17.6. For these fits we have altered the 
assumed stellar profiles from those used in Figure 8 to ones 
that produce better fits for low values of N ul . This 
approach amounts to decreasing the fluxes and broadening 
the profiles a bit. Despite these changes, Figure 10 shows 
that single-component fits with log N Hl = 17.6 are incapa- 
ble of fitting the data. The sides of the observed absorption 
lines are simply not steep enough. 

It seems sensible that hydrogen and deuterium should be 
well mixed and therefore share the same flow velocity and 
have line broadening parameters consistent with the same 
values of T and We therefore fit the H i and D i absorp- 
tion lines toward each star while imposing the constraint 
that the velocities and broadening parameters of H and D 
are consistent with the D i fits in Figure 6. The results are 
shown in Figure 11. The fits are again very poor. The H i 
absorption of the models are clearly blueshifted relative to 
the data, and the saturated cores of the model absorption 
profiles are clearly too narrow. Thus, for single-component 
models, the velocity and temperature of hydrogen cannot be 
forced to be consistent with the velocity and temperature of 
deuterium. 

3.2. T wo-Component Fits to the Line Profiles 

A second H i absorption component is obviously 
required to make the velocity and temperature of H i con- 
sistent with D I, Fe ii, and Mg n. This second component 
must have a H i column density high enough to modify the 
observed H i absorption significantly, but not high enough 
to produce observable D i, Fe ii, or Mg ii absorption. We 
first considered whether geocoronal H i absorption could 
be responsible. At the time of these observations, the Earth 
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Fig. 10.— Fits to the H I absorption lines of Cen A and a Cen B, in 
which the H I column density has been forced to be 10 17 6 . The assumed 
stellar profiles ( thin solid lines ) have been altered from those used in Fig. 8 
to ones that maximize the quality of fits with log N H1 = 17.6, but the fits 
( thick solid lines ) are very poor despite this change. This result proves that 
for single-component fits, log /V H] cannot be 17.6 and thus, D/H cannot be 
consistent with the D/H-value obtained for the Capella and Procyon lines 
of sight (Paper II). 

was moving toward a Cen A at 4- 7.3 km s“ 1 and toward a 
Cen B at +5.9 km s -1 in the heliocentric reference frame. 
Absorption by H i at these velocities could in principle 
produce the observed H i redshift, but these velocities are 
still well within the saturated core of the H i absorption. We 
find that for temperatures reasonable for geocoronal H i 
(T ~ 1000 K), H i column densities of about 10 17 cm 2 are 
necessary to broaden the absorption in the manner that we 
require. This column density is much too high, as typical H i 
column densities in the geocorona are 10 13 cm 2 . 
(According to the thermal exosphere model of Heiden 
(1988), typical values at 400 km are T = 1035 K and n H , = 
1 x 10 5 cm -3 .) 

Because the H i Lyman-a line is on the flat portion of the 
curve of growth, the addition of a small amount of hot H i 
in the line of sight will broaden the line core, and displacing 



1215.0 1215.2 1215.4 1215.6 1215.6 1216.0 1216.2 
Wavelength 

Fig. 1 1. — Fits to the H i and D i absorption lines of a Cen A and a Cen 
B in which the velocities and Doppler parameters of the lines are forced to 
be consistent with the velocities, temperatures, and nonthermal velocities 
measured from the D I fits in Fig. 6. The poor quality of the fits demon- 
strates that for a single-component fit, the velocity and temperature of H i 
cannot be consistent with those values for D I. 

this hot hydrogen component by a few km s -1 from the 
main component could explain the net redshift of the H i 
line relative to the D i line and the larger than expected 
value of b Hl . Thus, when we fit the D i and H i absorption 
lines simultaneously with two components and force consis- 
tency between D i and H I, the second component invari- 
ably ends up being hot and redshifted. These fits are shown 
in Figure 12, and the fitted parameters are given in Table 3. 
The assumed stellar profiles have been altered slightly from 
those used in Figure 8 to maximize the quality of the two- 
component fits. Since it is difficult to measure accurately the 
parameters of highly blended components such as those in 
Figure 12, we forced the velocities and Doppler parameters 
of the components representing the bulk of the ISM absorp- 
tion (dotted lines in the figure) to be consistent with the D i 
fits in Figure 6. Therefore, no errors are given for these 
parameters in Table 3. For the fits in Figure 12, we forced 
the D/H ratio to be the same for both components, but D i 


TABLE 3 


Interstellar Parameters for the Two-Component Fits to Lyman-oc 


Star 

Profile 

Model 

Velocity 
(km s -t ) 

^H! 

(km s *) 

log Nhi 

D/H 

(10 5 ) 

T 

x 2 

a Cen A 

1A 

-18.2 

9.3 

18.014 + 0.005 

0.59 + 0.03 

83766 

1.002 



-14.4 + 0.2 

20.01 + 0.67 

15.088 + 0.101 

0.59 + 0.03 

46.31 

1.002 

a Cen B 

IB 

-18.2 

9.7 

17.991 +0.003 

0.62 + 0.02 

75841 

1.589 



-13.8 + 0.1 

22.17 + 0.27 

14.860 ± 0.028 

0.62 + 0.02 

24.77 

1.589 

a Cen A 

2A 

-18.2 

9.3 

17.658 ± 0.014 

1.28 + 0.10 

36887 

1.105 



-15.7 + 0.2 

24.86 + 0.46 

14.843 + 0.034 

1.28 + 0.10 

21.21 

1.105 

a Cen B 

2B 

-18.2 

9.7 

17.597 + 0.007 

1.54 + 0.05 

30631 

1.436 



-15.3 + 0.1 

23.03 + 0.18 

15.007 + 0.017 

1.54 ±0.05 

33.39 

1.436 
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Fig. 12. — Two-component fits to the H i and D i absorption lines of a 
Cen A and a Cen B. The parameters of these fits are given in Table 3. The 
assumed stellar Lyman-a profiles ( thin solid lines) are nearly the same as 
those used in Fig. 8, although they have been altered slightly to improve 
the quality of the fits. The dotted lines show the component representing 
the bulk of the interstellar medium. For both the a Cen A and a Cen B fits, 
the velocity and Doppler parameter of this component has been forced to 
be consistent with the velocity and temperature of the ISM implied by the 
D i fits in Fig. 6. The dashed lines show a hot component that we believe 
represents absorption from hot H i near the heliopause (see § 4.5). The 
thick solid lines show the combination of the two absorption components, 
and the residuals of the fits are shown below each fit. 


absorption is completely negligible for the hot component 
(dashed lines in the figure) because the hot H i column 
density is so low. According to the fits, the second com- 
ponent has a Doppler parameter characteristic of tem- 
peratures of about 27,000 K and is redshifted by about 4 km 
s -1 relative to the main component. We believe that this 
second component is due to interstellar H i that has flowed 
into the heliosphere and is compressed and heated near the 
heliopause by interaction with the solar wind. The column 
density, redshift, and temperature of our second component 
are consistent with recent models of the outer heliosphere, 
which predict the existence of this “ hydrogen wall.” These 
models will be discussed further in § 4.5. 

We now consider whether the addition of a second com- 
ponent changes our conclusions about the D/H ratio for the 
a Cen line of sight. The D/H-values implied by the two- 
component fits in Figure 12 are not much different from the 
results of the single-component fits. This similarity is not 
surprising, considering that the stellar profiles are nearly the 
same as those in Figure 8. The question is this: Can we find 
an acceptable fit to the data with a completely different 
D/H-value by altering the assumed profiles? The answer 
was “no” for single-component fits (see Fig. 10). However, 
in Figure 13 we use the same stellar profiles that were used 
in Figure 10, and we find that acceptable two-component 



1215.0 1215.2 12154 1215.6 1215.8 1216.0 12162 

Wavelength 

Fig. 13. — Same as Fig. 12, but using different stellar emission profiles. 
For these fits, the derived interstellar H i column density is much lower 
than for the fits in Fig. 12 (see Table 3). 

fits do exist when using those profiles. On the basis of the 
fit parameters given in Table 3, the D/H-values for the 
fits in Figure 13 are much higher than the D/H — 
(5.7 + 0.2) x 10~ 6 value we obtained from the single- 
component fits. We have no direct way of determining 
whether the intrinsic profiles in Figure 12 (hereafter profiles 
1A and IB) or Figure 13 (profiles 2A and 2B) are better 
representatives of the stellar Lyman-a emission lines, and so 
we must consider both as plausible options. Thus, we find 
that the existence of the second component destroys our 
ability to measure precise values of D/H from our data. 

3.3. Three-Component Fits to the Line Profiles 

The two-component fits in Figures 12 and 13 are certain- 
ly improvements over the single-component fits in Figure 8, 
as demonstrated by both the residuals and the / 2 -values. 
Nevertheless, there is at least one remaining systematic dis- 
crepancy between the models and the data. All the fits in 
Figures 12 and 13 show a narrow flux excess near 1215.42 
A, which could be explained by temperature gradient in the 
ISM or temperature gradients in the heliospheric interface 
region (HIR), but we have represented the ISM and hydro- 
gen wall with single-temperature components. Another 
problem with the two-component models is that the a Cen 
line of sight passes not only through the Sun’s heliosphere, 
but also through an analogous “astrosphere” of a Cen. 
(Since the prefix “ helio-” implies a connection with the Sun, 
we replace the “helio-” prefix with “astro-.”) In other 
words, there should be two hydrogen walls along this line of 
sight, and it is important to determine how the existence of 
a second hydrogen wall may affect our analysis. 

According to Lallement & Bertin (1992), the interstellar 
flow vector for the solar neighborhood is directed toward 
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£ = 1 86?1 and b — — 16?4 with a velocity of v 0 = 25.7 km 
s" 1 . The a Cen line of sight is at an angle of 0 = 52° with 
respect to the upwind direction of this vector. In order to 
determine what the interstellar flow looks like in the rest 
frame of a Cen, we assume that the V G vector of Lallement 
& Bertin (1992) is appropriate (see also Lallement et al. 
1995), and we assume a distance, proper motion, and radial 
velocity for a Cen of d = 1.34 pc, = — 3'.'640 yr -1 , 
\i b = + 0'.'700 yr -1 , and v = —22.0 km s -1 (Hirshfeld, 
Sinnott, & Ochsenbein 1991). On the basis of these quanti- 
ties and a Cen’s position of a(2000) = 14 h 39 m 36 s and 
5(2000) = — 60°50'10", we find that the ISM flow vector for 
a Cen is (/ = 286?8, b = —78.7, i? 0 = 22.6 km s -1 ) and the 
line of sight toward the Sun is at an angle of 0 = 79° with 
respect to the upwind direction. 

The a Cen values for v 0 and 0 (22.6 km s -1 , 79 ) are not 
too different from the solar values (25.7 km s _1 , 52°). We 
also note that the stellar winds of a Cen A and a Cen B are 
probably similar to the solar wind, as a Cen A and a Cen B 
are both relatively inactive, main-sequence, solar-like stars. 
Thus, we conclude that the properties of a Cen’s astro- 
spheric interface region (AIR) along the Sun-a Cen line of 
sight are probably not that different from the properties of 
the solar HIR for this line of sight. Note that a Cen A and 
a Cen B probably share the same astrosphere, since the 
stars are only about 46 AU apart (Heintz 1982), and the 
solar example suggests distances of about 200 AU from the 
stars to the 44 astropause ” (the stellar equivalent of the 
heliopause). 

The ISM flow velocity for the a Cen line of sight is — 18.0 
km s“\ and the stellar radial velocity for a Cen is —22.0 
km s -1 . Models of the solar hydrogen wall (see, e.g., 
Baranov & Malama 1995) suggest that the velocity of a 
Cen’s hydrogen wall should be somewhere between —18.0 
km s _1 and —22.0 km s” 1 because the material in the wall 
should be decelerated relative to the star. Thus, the a Cen 
AIR absorption should produce either a slight blueshift of 
the H i absorption line relative to the D i, Mg n, and Fe n 
lines, or no observable shift, since the velocity separation 
between the ISM absorption and the a Cen AIR absorption 
is small. However, we observe a significant redshift for H i, 
not a small blueshift. In contrast, the solar HIR absorption 
should lie somewhere between 0 and — 18 km s“ 1 and can 
therefore produce a H I absorption line that is significantly 
redshifted relative to D i, Mg n, and Fe n, in agreement with 
our observations. Therefore, the hot components of Figures 
12 and 13 should definitely be attributed to the HIR rather 


than the AIR of a Cen. However, the fit parameters for this 
second component may be inaccurate because we have 
neglected the possible contributions of a Cen’s AIR to the 
H i absorption. 

We have performed three-component fits to the data to 
see how the addition of the third component (representing a 
Cen’s AIR) affects the parameters derived for the second 
component (representing the solar HIR). The results using 
profile models 1A and IB are shown in Figure 14, and 
Figure 15 shows the fits obtained assuming profiles 2A and 
2B. The parameters for all these fits are given in Table 4. It 
is very difficult, if not impossible, to separate absorption 
components as highly blended as the three components in 
Figures 14 and 15. To reduce the number of free parameters, 
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Fig. 14.— Same as Fig. 12, but with the addition of a third component 
(dash-dotted lines ) that is meant to represent absorption by hot hydrogen in 
the neighborhood of a Cen, analogous to the absorption near the Sun 
represented by the second component. The parameters of these fits are 
given in Table 4. 


TABLE 4 

Interstellar Parameters for the Three-Component Fits to Lyman-g 


Star 

Profile 

Model 

Velocity 
(km s“ ! ) 

(km s *) 

log N„, 

D/H 

(10 ~ s ) 

T 

X 2 

a Cen A 

1A 

-18.2 

9.3 

18.016 + 0.004 

0.58 + 0.03 

84131 

0.965 



-8.0 

19.30 + 0.60 

14.633 + 0.070 

0.58 + 0.03 

16.83 

0.965 



-20.0 

15.01 + 1.14 

15.555 ± 0.335 

0.58 ± 0.03 

180.6 

0.965 

a Cen B 

IB 

-18.2 

9.7 

17.997 ± 0.002 

0.60 + 0.02 

76816 

1.389 



-8.0 

21.35 + 0.23 

14.513 + 0.020 

0.60 + 0.02 

11.55 

1.389 



-20.0 

16.06 + 0.48 

15.171 +0.101 

0.60 + 0.02 

67.46 

1.389 

a Cen A 

2A 

-18.2 

9.3 

17.654 + 0.014 

1.29 + 0.08 

36549 

1.096 



-8.0 

22.39 + 0.56 

14.439 + 0.045 

1.29 ± 0.08 

9.29 

1.096 



-20.0 

24.83 ± 0.55 

14.575 + 0.034 

1.29 + 0.08 

11.47 

1.096 

a Cen B 

2B 

-18.2 

9.7 

17.603 ± 0.007 

1.52 + 0.06 

31028 

1.406 



-8.0 

22.00 + 0.27 

14.542 ± 0.020 

1.52 + 0.06 

11.98 

1.406 



-20.0 

21.29 + 0.26 

14.831 + 0.026 

1.52 ± 0.06 

24.09 

1.406 
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Fig. 15. - Same as Fig. 14, but using different stellar emission profiles. 
The assumed Lyman-a profiles are the same as those used in Fig. 13. 


we forced the velocity and Doppler parameter of the ISM 
component (dotted lines) to be consistent with the param- 
eters of the D i fits in Figure 6, as we did for the two- 
component fits. We also fixed the velocities of the HIR and 
AIR components to be -8 km s' 1 and —20 km s“\ 
respectively. These velocities are consistent with the theo- 
retical predictions of Baranov & Malama (1995). 

The derived parameters of the third component, rep- 
resenting a Cen’s AIR, are log = 15.03 ±0.43 and 
b\i\ = 19.5 + 4.6 km s \ respectively. These parameters 
varied quite a bit among the four fits listed in Table 4, which 
explains the large errors. The column density and Doppler 
parameter of the third component are similar to the values 
we find for the second component of the two- and three- 
component fits, consistent with our prediction that a Cen’s 
AIR should be similar to the HIR of the Sun. 

For the two-component fits, we found that the param- 
eters of the hot component representing the HIR are 
log 14.95 ±0.11, btf\ = 22.5 ±2.0 km s" 1 , and 

i’hi = — 14.8 ± 1.2 km s~ l . We now compare these param- 
eters with the parameters of the second component of the 
three-component fits in order to obtain our best estimates 
for the solar hydrogen wall parameters. Because we found it 
necessary to fix the velocities of the HIR and AIR com- 
ponents in Figures 14 and 15, we can only give a lower limit 
of > — 16 km s 1 for the velocity of the solar hydrogen 
wall. The addition of the third component slightly lowered 
the derived Doppler parameter of the HIR component from 
b[f i = 22.5 ± 2.0 km s~ l to 6 ( H 2 J = 21.3 ± 1.4 km s“\ and 
lowered the column density from log Ntf J = 14.95 ± 0.1 1 to 
log N\i \ = 14.53 ± 0.08. By averaging the results of all 8 
two- and three-component fits, we find that our best esti- 
mates and 1 o errors for the parameters of the HIR are 


i>h } > — 16 km s" 1 , btf 1 = 21.9 ±1.7 km s 1 
(corresponding to a temperature of T = 29,000 ± 5000 K), 
and log = 14.74 ± 0.24. 

4. DISCUSSION 

4.1. The Hi Column Density to a Cen 

We have found that the H i column density toward a Cen 
is determined mainly by fitting the shape of the observed 
Lyman-a line wings and thus depends critically upon the 
width of the wings of the assumed stellar emission line. On 
the other hand, the shape of the observed line cores pro- 
vides information on the hot gas that we think is near the 
heliopause (see § 4.5) rather than in the interstellar medium. 
Thus, the values of log /V HI and D/H depend upon whether 
profiles 1A and IB or profiles 2A and 2B are better repre- 
sentations of the intrinsic stellar emission lines. 

The wings of profiles 1A and IB are based on the mean 
quiet-Sun profile obtained with the HRTS (High Resolution 
Telescope Spectrograph) experiment and kindly provided 
to us by P. Brekke (Brekke et al. 1991). To our knowledge 
there are no high-resolution Lyman-a profiles from the 
whole Sun viewed as a star. The wings of profiles 2 A and 2B 
are roughly 10% broader than those of profiles 1A and IB. 
Profiles l A and IB represent our “best guesses” for the 
intrinsic stellar profiles, but are the broader 2A and 2B 
profiles sensible alternatives? HRTS spectra of the Sun with 
high spatial resolution show relatively narrow profiles over 
most of the solar disk but much brighter profiles with 
broader wings over active regions (see, e.g., Basri et al. 
1979). The fractional coverage of active regions on the Sun 
varies with the solar cycle, and the coverage of presumably 
analogous active regions on a Cen A and a Cen B is 
unknown. Active regions will broaden the Lyman-a profiles 
somewhat in spatially averaged stellar spectra. We expect 
that the Lyman-a profile for the Sun viewed as a star could 
have somewhat broader wings than the mean quiet-Sun 
profile used to construct profiles 1A and IB and that the 
same should be true for the a Cen stars. Thus, profiles 2A 
and 2B are plausible representations of the intrinsic emis- 
sion lines, and the interstellar hydrogen column density 
could be near log N Hl = 17.6 rather than 18.0 derived by 
using profiles 1A and IB. Because the error in the hydrogen 
column is dominated by the uncertain width of the intrinsic 
emission line profile rather than the random errors in the 
data or the formal errors derived by the fitting program, we 
adopt log N hi — 17.8 ± 0.3 as a rough but generous alloca- 
tion for this systematic error that includes both possible 
solutions. 

Landsman et al. (1984, 1986) analyzed WE Lyman-a 
spectra and reanalyzed archival Copernicus spectra of 
the a Cen stars. For a Cen A they obtained log N Hl = 
17.08-17.92, and for a Cen B they found log N Hl < 17.78. 
These ranges for the hydrogen column are considerably 
larger than our findings, but they overlap the values we 
derive from GHRS spectra. Our derived range for the 
neutral-hydrogen column toward a Cen is also consistent 
with the value N H1 = 9 x 10 17 cm -2 that Mewe et al. (1995) 
found by modeling the Extreme Ultraviolet Explorer 
(EUVE) spectrum of a Cen. R. Mewe (1995, private 
communication) has stated that the uncertainty in their 
hydrogen column is N H , = (9 ± 1) x 10 17 cm 2 , which cor- 
responds to log N H , = 17.95 ± 0.05. Finally, Quemerais el 
al. ( 1 994) derived the H i density, n H , = 0.1 65 ± 0.035 cm 3 , 
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for the LISM gas flowing into the solar system from com- 
parisons of radiative-transfer calculations with interplan- 
etary and geocoronal diffuse Lyman-a emission. This result 
corresponds to a column density of log N Hl = 17.83 ± 0.09 
toward a Cen, consistent with our results. Our value of 
log iV H i = 17.80 ± 0.30 corresponds to n H j = 0.076 0.30 
cm" 3 and D/H = (0.5 1.9) x 10” 5 . Our parameters are 
compared in Table 5 with those derived for the lines of sight 
to Procyon and Capella (Paper II). 

4.2. The Flow Velocity Along the a Cen Line of Sight 

We have found that the main velocity component along 
the line of sight toward a Cen (obtained from the D i, Mg n, 
and Fe n lines) has a bulk heliocentric flow speed of 
v = — 18.0 + 0.2 km s“ which is consistent with the speed 
of — 18.1 km s” 1 predicted by the flow vector in the Galac- 
tic center direction ( V G ) proposed by Lallement & Bertin 
(1992). If we had used their anti-Galactic center vector 
(F AG ), the projected velocity would have been —15.7 km 
s -1 . For Capella, Procyon, and many other stars, ISM 
absorption consistent with the F AG vector is observed, but 
no absorption is seen that is at the velocity predicted by the 
V G vector (see Papers I and II; Lallement et al. 1995). 
However, most or all of the line of sight toward a Cen 
appears to be through the V G gas. 

Which vector represents the true LISM flow vector for 
the solar neighborhood? Lallement et al. (1995) have 
argued that the Sun is located in gas that is flowing with the 
F AG vector, partly because a Cen seems to be the exception 
rather than the rule and partly because the 26 ± 1 km s" 1 
velocity of interstellar He i atoms measured by Ulysses 
(Witte et al. 1993) agrees with the modulus of the F AG vector 
rather than the faster velocity of 29 km s 1 predicted by the 
V G vector. These arguments do not explain, however, why 
we see no evidence of any material toward a Cen moving at 
the velocity predicted by the F AG vector. Assuming that the 
Sun does in fact lie in the AG cloud, Lallement et al. (1995) 
have estimated that the AG cloud can at most account for 
only about 2%-4% of the interstellar column density 
observed toward a Cen, and they have interpreted this small 
percentage to mean that the boundary between the AG and 
the G clouds (assuming they are in fact separate clouds) 
must be very nearby along the a Cen line of sight. 


4.3. The Temperature and Nonthermal Gas Motions 
Toward a Cen 

The mean temperature and turbulent velocity of the gas 
toward a Cen are T = 5400 + 500 K and £ = 1.20 ± 0.25 
km s" 1 , respectively. The nonthermal motions, which we 
have modeled as Gaussian microturbulence, are the same as 
for the Procyon (f = 1.21 ± 0.27 kms" 1 ) line of sight and 
are consistent with the value = 1.6 + 0.4 km s " 1 ) derived 
for the Capella line of sight (see Paper II). We have no way 
of knowing whether these nonthermal motions are a conse- 
quence of unresolved discrete velocity components, rotation 
of the LIC, expansion of the cloud along the line of sight, or 
true turbulent motions. The motions are definitely subsonic, 
and we conclude that there is no shock front in the LIC 
near the Sun and thus no shock heating at this time. 

On the other hand, the mean temperature of the inter- 
stellar gas toward a Cen is apparently lower than toward 
Capella (T = 7000 ± 500 K) and Procyon (T = 6900 ± 300 
K), both of which lie either outside or near the edge of the 
LIC. Lallement et al. (1995) also noted that the interstellar 
Fe n and Mg n lines toward a Cen A are narrow, which 
suggests that the gas in this line of sight is probably cooler 
than the gas in the LIC, but they did not have profiles of the 
D i line to determine the temperature accurately. Our data 
provide the first conclusive evidence that the interior of the 
G cloud is cooler than the edge of the LIC, which can be 
heated more effectively by the external UV radiation field 
and thermal conduction from the presumably hotter sur- 
rounding gas. 

The H i number density (n H1 = 0.076-0.30 cm" 2 ) and 
temperature (T = 5400 ± 500 K) of the interstellar gas 
toward a Cen correspond to a thermal pressure (including 
10% helium) P/k = 410-1950 cm" 3 K. We assume that this 
pressure is representative of the gas located toward the 
center of the G cloud. On the other hand, the Procyon and 
Capella data indicate a slightly higher temperature 
(T = 7000 K) when looking through the LIC in the 
outward direction. The density in the outward direction is 
not known, but if we assume n HI = 0.076-0.30 cm 2 , the 
thermal pressure P/k = 580-23 10 cm” 3 K. 

Wolfire et al. (1995) have computed models of the ISM in 
ionization and thermal equilibrium. They have found that 
three stable phases are possible for a limited range of pres- 


TABLE 5 


Comparison of Interstellar Parameters 


Parameter 

a Cen A/B 

(Profile 1) (Profile 2) 

Capella 

(Linsky et al. 1995) 

Procyon 

(Linsky et al. 1995) 

l 

316° 

316° 

163° 

214° 

b 

-1° 

-1° 

5° 

13° 

d (pc) 

1.34 

1.34 

12.5 

3.5 

v m . iR (km s" 1 ) 

-18.0 

-18.0 

22.0 + 0.9 

20.8 + 1.5 

^LISM S *) 

— 18.1 

-18.1 

22.0 

19.8 

"Wr (km s" 1 ) 

— 14 1 

-15.5 


23.4 + 1.5 

T (K) 

5400 + 500 

5400 + 500 

7000 + 500 

6900 ± 300 

£ (km s l ) 

1.20 + 0.25 

1.20 + 0.25 

1.6 + 0.4 

1.21 + 0.27 

N„, (10 17 cm' 2 ) 

10.1 

4.2 

17.1-17.6 

11.2-11.8 

<« HI )(cm' 3 ) 

0.24 

0.102 

0.0443-0.0456 

0.104-0.109 

(PC) 

D/H (10" 5 ) 



1.9— 7.3 

1. 2-4.8 

0.61 

1.41 

1.60 + 0.09 

1.6 

log (Mg n/H i) 

-5.28 

-4.91 

-5.41 

-5.50 

log (Fe ii/H i) 

-5.55 

-5.18 

-5.75 

-5.76 

D i/Mg ii 

1.2 + 0.2 

1.2 + 0.2 

4.1 + 0.4 

5.1 +0.7 
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sure: (1) a hot ionized medium (HIM) with T ~ 10 6 K, (2) a 
warm neutral medium (WNM) or warm ionized medium 
(WIM) with T ~ 8000 K, and (3) a cold neutral medium 
(CNM) with T ^ 50 K. The LIC is presumably in the 
WNM phase, in which heating is dominated by photoelec- 
tric heating of polycyclic aromatic hydrocarbons (PAHs) 
and cooling is dominated by radiation in the Lyman-a, C n 
(158 //m), and O i (63 pm) lines, and recombination to posi- 
tive ions. If P/k < 990 cm” 3 K, the HIM and WNM are 
stable, but the CNM is not stable, and there is no evidence 
for the CNM phase in the LIC. For their assumed standard 
conditions of the incident far ultraviolet radiation field and 
abundances and for /V H[ = 10 18 cm” 2 , their equilibrium 
model has T = 9200 K when n H , = 0.15 cm ” 3 . For equi- 
librium conditions and n Hl = 0. 1 0 cm " 3 , significantly lower 
pressures and temperatures are possible for lower values of 
the dust/gas (D/G) ratio; P/k = 580 cm ” 3 K corresponds to 
D/G ^0.2 of the nominal values cited in their paper. WNM 
temperatures are very insensitive to abundances and the 
FUV radiation field when n Hl = 0.10 cm 3 . Alternatively, 
the LIC may not be in ionization and thermal equilibrium. 
At present we cannot say which is the correct explanation 
for the low temperature and pressure in the LIC. 

4.4. Gas-Phase Abundances o/ Fe and Mg for the 
ol Cen Line of Sight 

We now estimate the abundances of Fe and Mg in the a 
Cen line of sight by making the same assumptions as in 
Papers I and II. For the column densities of Fe n and Mg n 
we average the values for a Cen A and a Cen B: log N hc n = 
12.45 and log N MgII = 12.72. For log N HI we have adopted 
17.80 ± 0.30. Therefore, log (N hel] /N Hl ) = -5.05 to -5.65 
and log (/V MgII /N HI ) = -4.78 to -5.38. 

Comparing these numbers to the solar abundances of Fe 
and Mg, log (Fe/H) 0 = —4.49 and log (Mg/H) 0 = —4.41 
(Anders & Grevesse 1989), and assuming that the contribu- 
tions to Fe ii and Mg n from ionized gas are negligible, we 
derive the logarithmic depletions, Z)(Fe) = (N Fe „/N H1 ) — 
log (Fe/H) G = -0.56 to -1.16 and Z)(Mg) = (JV MgI1 //V H ,) 
— log(Mg/H) 0 = —0.37 to —0.97 in the warm gas toward 
a Cen. It appears that Fe and Mg are less depleted toward a 
Cen than they are toward Capella [D(Fe) = —1.28 and 
D(Mg) = —1.03; Paper I] and Procyon [D(Fe) = —1.27 
and D( Mg) = — 1.09; Paper II]. 

The a Cen Mg and Fe depletion values can be compared 
to Copernicus estimates of D w (Mg) = —0.3 and D w (Fe) = 
— 1.4 in the warm interstellar medium (Jenkins, Savage, & 
Spitzer 1986). The depletion for Fe is smaller than our a Cen 
measurement, and the value for Mg is somewhat larger. 
Most of the difference for Mg may be due to a systematic 
error in the Rvalues for the far-ultraviolet Mg n lines near 
1240 A used in the Copernicus study. Sofia, Cardelli, & 
Savage (1994) found that the far-UV lines of Mg n yield 
Mg ii column densities 0.67 dex larger than those obtained 
from an analysis of the very strong but damped Mg ii lines 
near 2800 A. Because the /-values of the Mg ii 2800 A lines 
are well determined, Sofia et al. (1994) recommended a 0.67 
dex revision to the /-values for the Mg n far-UV doublet. 
Although there may be difficulty with the Mg n result, the 
measurements for Fe n clearly reveal that substantial gas- 
phase depletion occurs in the cloud toward a Cen. Evidently 
the Fe in the cloud is found mostly in interstellar dust. 

Because of the difficulty in measuring the column density 
of H i, the D/H ratio and the depletions of Mg and Fe are 


highly uncertain. Therefore, we are unable to determine 
whether D/H varies in the LISM, but we do find that the 
Mg and Fe depletions are at least slightly low'er than those 
found toward Capella and Procyon. The D i/Mg ii ratios, 
which are independent of the uncertain hydrogen column, 
are shown in Table 5 for a Cen, Capella, and Procyon. We 
find that D i/Mg ii is nearly a factor of 4 smaller for the a 
Cen line of sight than for the Capella and Procyon lines of 
sight. Thus , our data conclusively demonstrate that ISM 
abundances do vary significantly over distance scales of only a 
few parsecs , although we are not certain whether it is D i or 
Mg Ii/Fe ii that is actually varying. Because the Mg ii and 
Fe n depletions in the interstellar medium arc known to be 
variable (Jenkins 1987), the Mg n and Fe ii ions seem to be 
the most likely culprits. This argument favors the use of 
Lyman-a profiles 2A and 2B rather than profiles 1 A and 1 B, 
because the fits with the former profiles suggest hydrogen 
columns consistent with variable Mg ii and Fe ii depletions 
and constant D/H in the LISM. 

4.5. A New Test of Gasdynamical Models of the Interstellar 
Gas and Solar Wind Interaction: Evidence for Compressed , 
Heated Interstellar Neutral Hydrogen 

For the past 25 years a succession of space instruments 
have studied the diffuse Lyman-a emission produced when 
H i atoms in the heliosphere (the region of space dominated 
by the solar wind) resonantly scatter solar Lyman-a 
photons. For example, Quemerais et al. (1994) inferred a 
value for the interplanetary H i density far from the Sun, 
n H1 = 0.165 + 0.035 cm” 3 , by analyzing Lyman-a absorp- 
tion cell count rates obtained with the Atmospheric Lyman 
Alpha Experiment flown on the Space Shuttle in 1992. 
Their paper also provides a summary of the earlier obser- 
vations of diffuse Lyman-a emission. 

High-resolution spectra of the diffuse Lyman-a emission 
obtained with the GHRS and earlier instruments provide 
two critical facts concerning the properties of the H i atoms 
in the heliosphere. First, the heliocentric inflow speed is 
20-21 km s” 1 (first measured by Adams & Frisch 1977) 
rather than the 26 km s” 1 velocity that characterizes both 
the LIC gas beyond the heliosphere and the neutral He 
inside the heliosphere (Lallement & Bertaux 1990; Lalle- 
ment, Bertaux, & Clarke 1993; Clarke et al. 1995). The 
direction of the H i flow in the heliosphere and in the LIC 
are nearly identical (Lallement & Bertin 1992). Thus, H i is 
decelerated somewhere in the heliosphere by 5 6 km s” 1 
relative to the flow in the LIC. Second, the profile of the 
diffuse emission in both the upwind (the direction from 
which the interstellar gas appears to flow) and downwind 
directions can be fitted by a 30,000 K Voigt profile rather 
than the ^7000 K temperature generally thought to char- 
acterize the interstellar gas, which implies that the H i is 
either heated or that the interplanetary flow pattern is 
complex (Clarke et al. 1995). 

These spectroscopic data have stimulated the computa- 
tion of self-consistent multifluid gasdynamical models of the 
interaction between the inflowing interstellar gas and the 
expanding solar wind. These models include charge- 
exchange collisions in which interstellar H i atoms exchange 
an electron with solar wind protons, forming secondary H i 
atoms. Without charge exchange, the primary H i atoms 
from the interstellar medium would penetrate unimpeded 
through the solar system, as previously thought, except for 
photoionization and radiation pressure within 30 AU. The 
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original models of Baranov & Malama (1993), which were 
for only one set of interstellar parameters (n H j = 0. 14 cm “ 3 , 
n p = 0.07 cm -3 ), have now been extended by Baranov & 
Malama (1995), Williams et al. (1996), and Zank et ah (1996) 
to include a range of interstellar parameters. 

These new gasdynamical models typically show three dis- 
continuities in the gas flow. The outermost discontinuity is 
a bow shock located at 200-500 AU from the Sun 
(depending on n p in the ISM) in the upwind direction, where 
the flow of the ionized (plasma) component of the inter- 
stellar gas changes direction. The innermost discontinuity is 
the termination shock located at 100-200 AU, where the 
solar wind changes from a supersonic to a subsonic flow. 
Between these two shocks is the heliopause, the contact 
surface that separates the flow of interstellar plasma from 
that of the solar wind plasma. The location of these discon- 
tinuities depends critically on the uncertain fractional ion- 
ization of the interstellar medium; higher ionization leads 
to the boundaries lying closer to the Sun. In earlier helio- 
spheric models the inflowing H i does not see these shocks, 
but rather flows through the expanding solar wind without 
collisions until about 30 AU from the Sun where photoion- 
ization and radiation pressure produce a cavity of neutral 
H i. In such models, beyond 30 AU from the Sun the density 
of H i is nearly constant and equal to the interstellar 
density, and the H i flow speed is also the same as the 
interstellar flow. 

The main difference between the new and old models is 
the inclusion of charge-exchange collisions between the 
inflowing H i and the expanding solar wind protons, which 
then become secondary H i atoms with different momenta 
than the primaries. In these models, interactions between 
the primary and secondary H i atoms compress, decelerate, 
and heat the inflowing H I near the heliopause. This pile-up 
region of compressed H i in the upwind direction (first 
shown by Baranov, Lebedev, & Malama 1991) is now called 
the “hydrogen wall.” In the models of Baranov & Malama 
(1995), the temperature of the hydrogen wall gas is about 
20,000 K, whereas the models of Williams et al. (1996) indi- 
cate that 35,000 K is a typical upwind temperature for this 
gas all of the way into the inner solar system. Typical 
column densities for this warm H i are 10 15 cm' 2 both for 
typical LIC pressures and when the pressure is artificially 
raised by a factor of 3.5 to eliminate the bow shock. This 
higher pressure simulates the uncertain extra pressure terms 
from interstellar magnetic fields and cosmic rays. In the 
hydrogen wall the flow decelerates from the interstellar 
value of 26 km s" 1 to 7-23 km s' 1 depending on the 
assumed degree of ionization of the interstellar gas — higher 
ionization leads to more deceleration and forces the wall 
closer to the Sun. These models also predict 10 6 K (hot) H i 
primarily in the downwind direction. Charge exchange is 
unimportant for He I, so the He i atoms throughout the 
heliosphere retain the density and flow vector of the LIC. 

Although recent models of the solar wind and interstellar 
medium interaction are rather sophisticated, more realistic 
extensions of this theoretical work are needed and would 
benefit from clear observational constraints on the impor- 
tant parameters. Quemerais et al. (1995) searched for an 
increase in the diffuse Lyman-a emission in the upwind 
direction that would be a signature of the hydrogen wall in 
their analysis of low-dispersion spectra from the Voyager 1 
and Voyager 2 spacecraft, which were then 54 and 40 AU 
from the Sun. They indeed found an increase in the diffuse 


Lyman-a emission in the upwind direction, but they were 
unable to conclusively determine whether this increase was 
due to an enhancement of the H i density predicted by the 
models or to Galactic diffuse Lyman-a emission. GHRS 
Lyman-a spectra for the line of sight toward Mars show 
broad diffuse emission at radial velocities expected for the 
decelerated H i flow in the inner solar system. Clarke et al. 
(1995) suggested that the broad width of this emission could 
be due to a dispersion in the flow velocities along this line of 
sight, but the broad widths could be simply interpreted as 
thermal broadening with T = 30,000 K, which would be 
consistent with the theoretical models (see, e.g., Williams et 
al. 1996). 

Our analysis of the Lyman-a absorption along the line of 
sight to a Cen led us to include a warm component with the 
parameters T = 29,000 ± 5000 K, log N Hl = 14.74 + 0.24, 
and a flow velocity of v H , > — 16 km s" *, which suggests at 
least a small deceleration relative to the LIC inflow velocity 
at an angle of 52° relative to the upwind direction (see § 3.3). 
These parameters are within the range of recent theoretical 
models and support the basic assumptions and results of the 
models. 

The importance of our measurement of the absorption by 
warm hydrogen toward a Cen is that for the first time we 
are measuring the column density of H i along the entire 
line of sight through the heliosphere, rather than an in situ 
measurement of scattered H i flux, which may not sample 
well the H i located in the heliosphere far from the Sun 
where the solar Lyman-a flux is low. Also, there is no confu- 
sion with possible diffuse Galactic emission or uncertainties 
concerning the solar Lyman-a flux. Our data provide the 
first direct measurement of the temperature of the com- 
pressed and heated interstellar H i in the hydrogen wall that 
can be used to constrain the theoretical models, a Cen is 
particularly well suited for studying the hydrogen wall 
because (1) the stars lie in the upwind direction where the 
wall is predicted to be (whereas stars previously observed by 
the GHRS all lie in the downwind direction); (2) the inter- 
stellar hydrogen column is very low, so the interstellar 
absorption is narrow enough to not obliterate absorption 
by warm hydrogen in the wall; and (3) a Cen A and a Cen B 
are bright enough in the UV to provide high S/N spectra. 
Lyman-a spectra for other lines of sight should be studied to 
confirm our conclusions. 

4.6. Astronephography of the LIC 

While a Cen is probably located inside the G cloud, the 
other stars observed by the GHRS, such as Procyon and 
Capella, probably lie outside the LIC, which is also a warm, 
neutral cloud (Lallement et al. 1995). We can use the know- 
ledge that a Cen lies within the G cloud to estimate the 
mean H i density along this line of sight. Our value of 
log Nhi = 17.80 + 0.30 corresponds to n Hl = 0.15 cm' 3 
with an uncertainty of a factor of 2. To our knowledge, this 
is the first direct measurement of the H i density within a 
local cloud, rather than an inference of n H , on the basis of 
Lyman-a scattered emission inside the heliosphere and 
theoretical models of the solar wind and interstellar gas 
interaction. 

If we adopt the value 0.15 cm" 3 as representative of the 
H i density in the G cloud and the LIC, we can estimate the 
distance to the edge of the LIC (d edge ) in the direction of 
Procyon, Capella, and other stars. We find that the distance 
to the edge of the LIC is only 3.9 pc toward Capella and 
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only 2.5 pc toward Procyon. Since the distance to Procyon 
is 3.5 pc, there is room for the second cloud detected along 
this line of sight (see Paper II) to lie outside of the LIC. 
These values for d edge are uncertain because of the factor of 
2 uncertainty in the H i density in the G cloud and the 
assumption that the densities of the G cloud and the LIC 
are the same. 

Accurate measurements of the H i column toward other 
stars obtained with the GHRS are now becoming available. 
For example, Gry et al. (1995) found that component 1 in 
the line of sight toward e CMa has the local cloud velocity 
and log A h , = 17.34. The hot white dwarf G191— B2B 
shows a component (called cloud B) with the local cloud 
velocity and log N Hl = 18.27 (Lemoine et al. 1995). Also, 
Frisch (1994) has estimated that component 1 along the line 
of sight to Sirius is mostly ionized with a total column of 
5.1 x 10 17 cm' 1 2 . By using the same technique as before, we 
estimate the distance to the cloud edge toward each star to 
be 0.47 pc (e CMa), 4.0 pc (G191 — B2B), and 1.1 pc (Sirius). 
Since Sirius is at a distance of 2.7 pc, there is room for the 
second cloud along its line of sight (see Lallement et al. 
1994) to lie outside of the LIC. 

With five distances to the edge of the LIC in different 
directions, we can begin to draw the surface of the local 
cloud. Frisch (1994) and Lallement et al. (1995) began this 
morphological study by using inferred rather than mea- 
sured values of n HI . As more diverse lines of sight are 
studied, we can construct a three-dimensional representa- 
tion of the cloud. Since this is a new field of study, it requires 
a new name. In analogy with “geography,” we propose the 
name “ astronephography ” based on the Greek word vecpoa 
(nephos), meaning “ cloud ” 

5. CONCLUSIONS 

We summarize here the results of our analysis of GHRS 
echelle spectra of interstellar absorption in the short (1.34 
pc) lines of sight toward a Cen A and a Cen B : 

1. The interstellar absorption features detected in the 
resonance lines of Mg ii, Fe ii, and D i all have the same 
heliocentric velocity (t? = —18.0 + 0.2 km s' 1 ), which is 
consistent with the flow vector for the G cloud proposed by 
Lallement & Bertin (1992). The temperature and non- 
thermal velocity inferred from these line profiles are 
T = 5400 + 500 K and f = 1.20 + 0.25 km s' 1 , respec- 
tively. The temperature is significantly lower than the 7000 
K value previously obtained for the Capella and Procyon 
lines of sight that extend through the outer edge of the local 
interstellar cloud (LIC). The turbulent velocity toward a 
Cen is the same as that found for the Procyon line of sight 
and is very subsonic. The lower temperature toward a Cen 
could result from the G cloud being more shielded from 
ionizing flux from hot stars such as e CMa located outside 
of the LIC or perhaps could indicate a lower dust/gas ratio 
in the G Cloud. 

2. Our single-component fits to the H i lines suggest 

log N h , = 18.03 ± 0.01 and D/H = (5.7 ± 0.2) x 10" 6 . We 

had hoped that the very short lines of sight toward a Cen A 

and <x Cen B would contain gas with only one set of physical 
parameters (e.g., flow velocity, temperature, and turbulent 
velocity), but we found that the H i absorption is broader 

than expected from the other lines with a Doppler param- 


eter (b„ |= 11.80 km s' 1 ), which implies a significantly 
warmer temperature of 8350 K. Furthermore, the velocity 
of the H i absorption (v — — 15.8 + 0.2 km s ~ 1 ) is redshifted 
by about 2.2 km s' 1 with respect to the velocity of the Fe ii. 
Mg ii, and D i lines. 

3. The most sensible way to resolve the discrepancy 
between H i and the other lines is to add a second com- 
ponent to the H i lines. This second component is hotter 
( T % 30,000 K), is redshifted relative to the primary com- 
ponent by 2- 4 km s“ \ and has a column density too low to 
be detected in the Fe n, Mg ii, and D i lines. We propose 
that the gas responsible for this second component is 
located near the heliopause and is the heated neutral hydro- 
gen from the interstellar medium that is compressed by the 
solar wind. This so-called “hydrogen wall” is predicted by 
recent multifluid gasdynamical models of the interstellar 
gas and solar wind interaction. Our data provide the first 
measurements of the temperature and column density of H i 
in the hydrogen wall. Our estimates of the parameters of the 
solar hydrogen wall are log Ntf J = 14.74 + 0.24, 6tf) = 
21.9 ± 1.7 km s' 1 (corresponding to a temperature of 
T = 29,000 ± 5000 K), and rtf! > - 16 km s' 

4. Unfortunately, the existence of this heated H I com- 
promises our ability to compute the H i column density of 
the interstellar medium accurately because, with slight 
alterations to our assumed stellar Lyman-a emission line 
profiles, acceptable two-component fits also exist with 
log N h [ ~ 1 7.6. We, therefore, quote large error bars for the 
hydrogen column density along the a Cen line of sight: 
log JV HI = 17.80 ± 0.30, n Hl = 0.076-0.30 cm ' 3 , and 
D/H = (0.5 1.9) x 10" 5 . This range in D/H is consistent 
with the value D/H = 1.6 x 10“ 5 previously derived for the 
Capella and Procyon lines of sight. With the large range of 
acceptable D/H-values, we cannot tell whether the D/H 
ratio varies or is constant in the local interstellar medium, 
but we do find that the D i/Mg ii ratio for the a Cen line of 
sight is about 4 times smaller than for the Capella and 
Procyon lines of sight. Therefore, either D/H or the Mg 
depletion varies significantly over distance scales of only a 
few parsecs. The known variations in the Mg ii depletions in 
the interstellar medium suggest that the Mg ii depletion is 
different for the a Cen line of sight, in which case D/H 
would be the same as for the Capella and Procyon lines of 
sight. 

5. Since a Cen lies inside the G cloud, we estimate the 
mean density for this line of sight to be n HI = 0.15 cm 3 
with an uncertainty of a factor of 2. This is the first direct 
measurement of n H1 in a nearby interstellar cloud, since the 
other stars observed by the GHRS likely lie outside the 
LIC. We can use the H i density in the G cloud to estimate 
the distances to the edge of the LIC along the lines of sight 
to Capella, Procyon, Sirius, G191 — B2B, and e CMa. 
These distances provide a rough outline for the LIC and the 
beginning of “ astronephography ” as a field of research. 


This work is supported by NASA Interagency Transfer 
S-56460-D to the National Institute of Standards and Tech- 
nology. We wish to thank V. B. Baranov, Priscilla Frisch, 
Rosine Lallement, and Lance Williams for sharing their 
expertise and for sending to us preprints of their recent 
papers, which were very useful. We also thank Blair Savage 
for his helpful comments on the manuscript. 



270 


LINSKY & WOOD 


REFERENCES 


Adams, T. F., & Frisch, P. C. 1977, ApJ, 212, 300 
Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Baranov, V. B., Lebedev, M. G., & Malama, Y. G. 1991, ApJ, 375, 347 
Baranov, V. B., & Malama, Y. G. 1993, J. Geophys. Res., 98, no. A9, 15157 
. 1995, J. Geophys. Res., 100, no. A8, 14755 

Basri, G. S., Linsky, J. L., Bartoe, J.-D. F., Brueckner, G., & Van Hoosier, 
M.E. 1979, ApJ, 230, 924 
Brandt, J. C., et al. 1994, PASP, 106, 890 

Brekke, P., Kjeldseth-Moe, O., Bartoe, J.-D. F., & Brueckner, G. E. 1991, 
ApJS, 75, 1337 

Bruston, P., Audouze, J., Vidal-Madjar, A., & Laurent, C. 1981, ApJ, 243, 
161 

Cerruti-Sola, M., Cheng, C.-C., & Pallavicini, R. 1992, A&A, 256, 185 
Clarke, J. T., Lallement, R , Bertaux, J.-L., & Quemerais, E. 1995, ApJ, 448, 
893 

Dupree, A. K., Baliunas, S. L., & Shipman, H. L. 1977, ApJ, 218, 361 
Ferlet, R., Lecavelier, A., Vidal-Madjar, A., Bertin, P., Deleuil, M., Gry, C., 
Lagrange-Henry, A.-M., & Lallement, R. 1995, A&A, 297, L5 
Frisch, P. C. 1994, Science, 265, 1423 
. 1995, Space Sci. Rev., 72, 499 

Gilliland, R. L. 1994, GHRS Instrument Science Report 063 (Baltimore: 
STScI) 

Gry, C., Lemonon, L., Vidal-Madjar, A., Lemoine, M., & Ferlet, R. 1995, 
A&A, 302, 497 

Heap, S. R., et al. 1995, PASP, 107, 871 

Heiden, A. E. 1988, J. Geophys. Res., 88, no. A12, 10170 

Heintz, W. D. 1982, Observatory, 102, 42 

Hirshfeld, A., Sinnott, R. W., & Ochsenbein, F. 1991, Sky Catalogue 
2000.0, Vol. 1 (Cambridge: Cambridge Univ. Press) 

Jenkins, E. B. 1987, in Interstellar Processes, ed. D. H. Hollenbach & H. A. 

Thronson, Jr. (Dordrecht: Reidel), 533 
Jenkins, E. B., Savage, B. D., & Spitzer, L. 1986, ApJ, 301, 355 
Lallement, R., & Bertaux, J.-L. 1990, A&A, 23 1 , L3 
Lallement, R., Bertaux, J.-L., & Clarke, J. 1993, Science, 260, 1095 
Lallement, R., & Bertin, P. 1992, A&A, 266, 479 

Lallement, R., Bertin, P., Ferlet, R., Vidal-Madjar, A., & Bertaux, J. L. 
1994, A&A, 286, 898 


Lallement, R., Ferlet, R., Lagrange, A. M., Lemoine, M., & Vidal-Madjar, 

A. 1995, A&A, 304, 461 

Landsman, W. B., Henry, R. C., Moos, H. W., & Linsky, J. L. 1984, ApJ, 
285, 801 

Landsman, W. B., Murthy, J., Henry, R. C., Moos, H. W., & Russell, J. L. 
1986, ApJ, 303, 791 

Landsman, W., & Simon, T. 1993, ApJ, 408, 302 

Lemoine, M., Vidal-Madjar, A., Ferlet, R., Bertin, P., Gry, C., & Lallement, 
R. 1995, in The Light Element Abundances, ed. P. Crane (Berlin: 
Springer), 233 

Linsky, J. L., et al. 1993, ApJ, 402, 694 [Paper I] 

Linsky, J. L., Diplas, A., Wood, B. E., Brown, A., Ayres, T. R., & Savage, 

B. D. 1995, ApJ, 451, 335 [Paper II] 

Mewe, R., Kaastra, J. S., Schrijver, C. J., van den Oord, G. H. J., & Alke- 
made, F. J. M. 1995, A&A, 296, 477 

Quemerais, E., Bertaux, J.-L., Sandel, B. R.. & Lallement, R. 1994, A&A, 
290, 941 

Quemerais, E„ Sandel, B. R., Lallement, R., & Bertaux, J.-L. 1995, A&A, 
299, 249 

Robinson, R. D., Blackwell, J., Feggans, K., Lindler, D., Norman, D., & 
Shore, S. N. 1992, A User’s Guide to the GHRS Software, Version 2.0 
(Greenbelt : Goddard Space Flight Center) 

Soderblom, D. R., Hulbert, S. J., Leitherer, C., & Sherbert, L. E. 1994, HST 
Goddard High Resolution Spectrograph Instrument Handbook, 
Version 5.0 (Baltimore: STScI) 

Sofia, U. J., Cardelli, J. A., & Savage, B. D. 1994, ApJ, 430, 650 
Vladilo, G., Molaro, P., Crivellari, L., Foing, B. H., Beckerman, J. E., & 
Genova, R. 1987, A&A, 185, 233 

Williams, L. L., Hall, D. T., Pauls, H. L., & Zank, G. P. 1996, in Proc. Solar 
Wind 8 Workshop, in press 
Wilson, O. C., & Bappu, M. K. V. 1957, ApJ, 125, 661 
Witte, M., Rosenbauer, H., Banaszkewicz, M., & Fahr, H. 1993, Adv. Space 
Res., 13, no. 6, 121 

Wolfire, M. G., Hollenbach, D., McKee, C., Tielens, A. G. G. M., & Bakes, 
E. L. 0.1995, ApJ, 443,152 

Zank, G. P., Pauls, H. L., Williams, L. L., & Hall, D. T. 1996, in Proc. Solar 
Wind 8 Workshop, in press 



